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Abstract The great flexibility of a view camera allows the 
acquisition of high quality images that would not be possible 
any other way. Bringing a given object into focus is however 
a long and tedious task, although the underlying optical laws 
are known. A fundamental parameter is the aperture of the 
lens entrance pupil because it directly affects the depth of 
field. The smaller the aperture, the larger the depth of field. 
However a too small aperture destroys the sharpness of the 
image because of diffraction on the pupil edges. Hence, the 
desired optimal configuration of the camera is such that the 
object has the sharpest image with the greatest possible lens 
aperture. In this paper, we show that when the object is a 
convex polyhedron, an elegant solution to this problem can 
be found. The problem takes the form of a constrained op- 
timization problem, for which theoretical and numerical re- 
sults are given. 

Keywords Large format photography • Computational 
photography • Scheimpflug principle 



1 Introduction 

Since the registration of Theodor Scheimpflug's patent in 
1904 (see [6J), and the book of Larmore in 1965 where a 
proof of the so-called Scheimpflug principle can be found 
(see in p. 171-173]), very Uttle has been written about the 
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Fig. 1 The Sinar e (a) and its metering back (b). 



mathematical concepts used in modern view cameras, until 
the development of the Sinar e in 1988 (see Figure [T}. A 
short description of this camera is given in [T p. 23]: 

The Sinar e features an integrated electronic computer, 
and in the studio ojfers a maximum of convenience and opti- 
mum computerized image setting. The user-friendly software 
guides the photographer through the shot without technical 
confusion. The photographer selects the perspective ( cam- 
era viewpoint) and the lens, and chooses the areas in the 
subject that are to be shown sharp with a probe. From these 
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scattered points the Sinar e calculates the optimum position 
of the plane of focus, the working aperture needed, and in- 
forms the photographer of the settings needed 

Sinar sold a few models of this camera and discontinued 
its development in the early nineties. Surprisingly, there has 
been very little published user feedback about the camera 
itself. However many authors started to study (in fact, re- 
discover) the underlying mathematics (see e.g. lIU and the 
references therein). The most crucial aspect is the consider- 
ation of depth of field and the mathematical aspects of this 
precise point are now well understood. When the geomet- 
rical configuration of the view camera is precisely known, 
then the depth of field region (the region of space where ob- 
jects have a sharp image) can be determined by using the 
laws of geometric optics. Unfortunately, these laws can only 
be used as a rule of thumb when operating by hand on a 
classical view camera. Moreover, the photographer is rather 
interested in the inverse problem: given an object which has 
to be rendered sharply, what is the optimal configuration of 
the view camera? A fundamental parameter of this config- 
uration is the aperture of the camera lens. Decreasing the 
lens aperture diameter increases the depth of field but also 
increases the diffraction of light by the lens entrance pupil. 
Since diffraction decreases the sharpness of the image, the 
optimal configuration should be such that the object fits the 
depth of field region with the greatest aperture. 

This paper presents the mathematical tools used in the 
software of a computer controlled view camera solving this 
problem. Thanks to the high precision machining of its com- 
ponents, and to the known optical parameters of the lens and 
digital sensor, a reliable mathematical model of the view 
camera has been developed. This model allows the acqui- 
sition of 3D coordinates of the object to be photographed, as 
explained in Section |2l In Section [3] we study the depth of 
field optimization problem from a theoretical and numerical 
point of view. We conclude and briefly describe the archi- 
tecture of the software in Section |4] 

2 Basic mathematical modeling 

2.1 Geometrical model of the View Camera 

We consider the robotised view camera depicted in Figure 
|2a) and its geometrical model in Figure |2b). We use a 
global Euclidean coordinate system {0,Xi,X2,X-},) attached 
to the camera's tripod. The front standard, symbolized by its 
frame with center L of global coordinates L = (Li ,L2,i3)^, 
can rotate along its tilt and swing axes with angles Oi and 
^i. Most camera lenses are in fact thick lenses and nodal 
points H' and H have to be considered (see [5 | p. 43-46). 
The rear nodal plane, which is parallel and rigidly fixed to 
the front standard, passes through the rear nodal point H' . 



Sensor Plane 




Fig. 2 Geometrical model (a) and robotised view camera (b). 

Since L and H' do not necessarily coincide, the translation 
between these two points is denoted t*^. The vector HH' is 
supposed to be orthogonal to the rear nodal plane. 

The rear standard is symbolized by its frame with center 
S, whose global coordinates are given by S — {Sl,S2^S{)^ . 
It can rotate along its tilt and swing axes with angles Qg and 
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(ps- The sensor plane is parallel and rigidly fixed to the rear 
standard. The eventual translation between S and the center 
of the sensor is denoted by t^. 

The rear standard center S can move in the three Xi,X2 
and directions but the front standard center L is fixed. 
The rotation matrices associated with the front and rear stan- 
dard alt-azimuth mounts are respectively given by R'^ = 
R(0L, and = R(05, (j)s) where 

(COS0 — sin(^sin0 — sin^cos0' 
COS0 -sine 
— sin^ cos ^ sin cos cos 

The intrinsic parameters of the camera (focal length /, po- 
sitions of the nodal points H, H\ translations t^, t'^, image 
sensor characteristics) are given by their respective manu- 
facturers data-sheets. The extrinsic parameters of the cam- 
era are S, L, the global coordinate vectors of S and L, and 
the four rotation angles 9s, 0s, Ol, (j)^. The precise knowl- 
edge of the extrinsic parameters is possible thanks to the 
computer-aided design model used for manufacturing the 
camera components. In addition, translations and rotations 
of the rear and front standards are controlled by stepper mo- 
tors whose positions can be precisely known. In the follow- 
ing, we will see that this precise geometrical model of the 
view camera allows one to solve various photographic prob- 
lems. The first problem is the determination of coordinates 
of selected points of the object to be photographed. 

In the sequel, for sake of simplicity, we will give all al- 
gebraic details of the computations for a thin lens, i.e. when 
the two nodal points H, H' coincide. In this case, the nodal 
planes are coincident in a so-called lens plane. We will also 
consider that t"" = (0,0,0)^ so that L is the optical center of 
the lens. Finally, we also consider that t'^ = (0,0,0)^ so that 
S coincides with the center of the sensor 



Lens Plane 




Fig. 3 Graphical constmction of the image A of an object point X in 
the optical coordinate system. Green rays and blue rays respectively lie 
in the [L,xj,,X2) and {L,Xi,x\) planes and F is the focal point. 



Sensor Plane 



Lens Plane 




2.2 Acquisition of object points coordinates 

Let us consider a point X with global coordinates given by 
X — {X\^X2,X-},)^ . The geometrical construction of the im- 
age A of X through the lens is depicted in Figure |3] We have 
considered a local optical coordinate system attached to the 
lens plane with origin L. The local coordinates of X are given 
by X = (R*^)^' (X — L) and the focal point F has local coor- 
dinates (0,0,/)^. Elementary geometrical optics (see |5 | p. 
35-42) allows one to conclude that if the local coordinates of 
A are given by a = {a\,a2,a^)^ , then a^, X3 and / are linked 
by the thin lens equation given in its Gaussian form by 



1 

03 



1 

-^3 



1 



Since A lies on the {XL) line, the other coordinates are ob- 
tained by straightforward computations and we have the con- 



Fig. 4 Image formation when considering a lens with a pupil. 



jugate formulas 



/ 



^a. 

/ + fl3 



(1) 

(2) 



Bringing an object into focus is one of the main tasks of a 
photographer but it can also be used to calculate the coordi- 
nates of an object point. It is important to remember that all 
light rays emanating from X converge to A but pass through 
a pupil (or diaphragm) assumed to be circular, as depicted in 
Figure |4] Since all rays lie within the oblique circular cone 
of vertex A and whose base is the pupil, the image of X on 
the sensor will be in focus only if the sensor plane passes 
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through A, otherwise its extended image will be a blur spot. 
By using the full aperture of the lens, the image will rapidly 
go out of focus if the sensor plane is not correctly placed, 
e.g. by translating S into the xj direction. This is why auto- 
focus systems on classical cameras only work at near full 
aperture: the distance to an object is better determined when 
the depth of field is minimal. 

The uncertainty on the position of S giving the best focus 
is related to the diameter of the so-called "circle of confu- 
sion", i.e. the maximum diameter of a blur spot that is in- 
distinguishable from a point. Hence, everything depends on 
the size of photosites on the sensor and on the precision of 
the focusing system (either manual or automatic). This un- 
certainty is acceptable and should be negligible compared to 
the uncertainty of intrinsic and extrinsic camera parameters. 

The previous analysis shows that the global coordinates 
of X can be computed, given the position (m, v)^ of its image 
A on the sensor plane. This idea has been already used on the 
Sinar e, where the acquisition of (m, v)^ was done by using 
a mechanical metering unit (see Figure [T](b)). In the system 
we have developed, a mouse click in the live video window 
of the sensor is enough to indicate these coordinates. Once 
(m, v)^ is known, the coordinates of A in the global coordi- 
nate system are given by 

and its coordinates in the optical system by 

a=(RL)-^(A-L). 

Then the local coordinate vector x of the reciprocal image 
is computed with (|2|i and the global coordinate vector X is 
obtained by 

X = L + R^x. 

By iteratively focusing on different parts of the object, the 
photographer can obtain a set of points = {X^ ,. .. ,X"}, 
with « > 3, which can be used to determine the best configu- 
ration of the view camera, i.e. the positions of front and rear 
standards and their two rotations, in order to satisfy focus 
requirements. 

3 Focus and depth of field optimization 

In classical digital single-lens reflex (DLSR) cameras, the 
sensor plane is always parallel to the lens plane and to the 
plane of focus. For example, bringing into focus a long and 
flat object which is not parallel to the sensor needs to de- 
crease the aperture of the lens in order to extend the depth 
of field. On the contrary, view cameras with tilts and swings 
(or DLSR with a tilt/shift lens) allow to skew away the plane 
of focus from the parallel in any direction. Hence, bringing 



into focus the same long and flat object with a view camera 
can be done at full aperture. This focusing process is unfor- 
tunately very tedious. However, if a geometric model of the 
camera and the object are available, the adequate rotations 
can be estimated precisely. In the next sections, we will ex- 
plain how to compute the rear standard position and the tilt 
and swing angles of both standards to solve two different 
problems: 

1. when the focus zone is roughly flat, and depth of field 
is not a critical issue, then the object plane is computed 
from the set of object points If n = 3 and the points 
are not aligned then this plane is uniquely defined. If 
n > 3 and at least 3 points are not aligned, we compute 
the best fitting plane minimizing the sum of squared or- 
thogonal distances to points of Then, we are able to 
bring this plane into sharp focus by acting on: 

(a) the angles 9i and 0^ of the front standard and the 
position of the rear standard, for arbitrary rotation 
angles 0s, (j)s. 

(b) the angles 9s, (ps and position of the rear standard, 
for arbitrary rotation angles 9l, (in this case there 
is a perspective distortion). 

2. when the focus zone is not flat, then the tridimensional 
shape of the object has to be taken into account. 

The computations in case [la] are detailed in Section 13.11 
In Section 13.31 a general algorithm is described that allows 
the computation of angles 9l and 0/, of the front standard 
and the position of the rear standard such that all the object 
points are in the depth of field region with a maximum aper- 
ture. We give a theoretical result showing that the determi- 
nation of the solution amounts to enumerate a finite number 
of configurations. 

3.1 Placement of the plane of sharp focus by using tilt and 
swing angles 

In this section we study the problem of computing the tilt 
and swing angles of front standard and the position of the 
rear standard for a given sharp focus plane. Although the 
underlying laws are well-known and are widely described 
(see ||4]|8]|2l), the detail of the computations is always done 
for the particular case where only the tilt angle 9 is consid- 
ered. Since we aim to consider the more general case where 
tilt and swing angles are used, we will describe the various 
objects (lines, planes) and the associated computations by 
using linear algebra tools. 

3.1.1 The Scheimpflug and the Hinge rules 

In order to explain the Scheimpflug rule, we will refer to 
the diagram depicted in Figure |5] The Plane of sharp focus 
(abbreviated SFP) is determined by a normal vector n^^ and 
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Sensor Plane (SP) Lens Plane (LP) 



Shaq) Focus Plane (SFP) 




Scheimpflug Line (SL) 



Fig. 5 Illustration of the Scheimpflug rule. 



Parallel to Sensor 

Lens Plane Front Focal Plane 
SP (PSLP) LP (FFP) 




Hinge Line (HL) 



Fig. 6 Illustration of the Hinge rule. 

a point Y. The position of the optical center L and a vector 
normal to the sensor plane (abbreviated SP) are known. 
The unknowns are the position of the sensor center S and a 
vector n'^ normal to the lens plane (abbreviated LP). 

The Scheimplug rule stipulates that if SFP is into focus, 
then SP, LP and SFP necessarily intersect on a common line 
called the "Scheimpflug Line" (abbreviated SL). The dia- 



gram of Figure|5^ should help the reader to see that this rule 
is not sufficient to uniquely determine n'" and SP, as this 
plane can be translated toward if is changed accord- 
ingly. 

The missing constraints are provided by the Hinge rule, 
which is illustrated in Figure|6] This rule considers two com- 
plimentary planes: the front focal plane (abbreviated FFP), 
which is parallel to LP and passes through the focal point 
F, and the parallel to sensor lens plane (abbreviated PSLP), 
which is parallel to SP and passes through the optical center 
L. The Hinge Rule stipulates that FFP, PSLP and SFP must 
intersect along a common line called the Hinge Line (abbre- 
viated HL). Since HL is uniquely determined as the intersec- 
tion of SFP and PSLP, this allows one to determine n^, or 
equivalently the tilt and swing angles, such that FFP passes 
through HL and F. Then SL is uniquely defined as the in- 
tersection of LP and SFP by the Scheimpflug rule (note that 
SL and HL are parallel by construction). Since is already 
known, any point belonging to SL is sufficient to uniquely 
define SP. Hence, the determination of tilt and swing angles 
and position of the rear standard can be summarized as fol- 
lows: 

1 . determination of HL, intersection of FFP and SFP, 

2. determination of tilt and swing angles such that HL be- 
longs to FFP, 

3. determination of SL, intersection of LP and SFP, 

4. translation of S such that SL belongs to SP. 

3.1.2 Algebraic details of the computations 

In this section the origin of the coordinate system is the op- 
tical center L and the inner product of two vectors X and Y 
is expressed by using the matrix notation X^Y. All planes 
are defined by a unit normal vector and a point in the plane 
as follows: 



SP 
PSLP: 
LP: 
FFP 
SFP: 



{xe 
{xe 
{xe: 
{xe: 
{xe: 



(x-s)V = o} 



X^n^ = 



X^n^ = o} , 
xTnL-/ = o}, 
(X-Y)"^n^'' = o}, 



where the equation of FFP takes this particular form because 
the distance between L and F is equal to the focal length 
/ and we have imposed that nk, > 0. The computations are 
detailed in the following algorithm: 

Algorithm 1 

Step 1 : compute the Hinge Line by considering its para- 
metric equation 

HL= {XeM^ 3teR, X = W + fV}, 
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where V is a direction vector and W is the coordinate 
vector of an arbitrary point of HL. Since this Une is the 
intersection of PSLP and SFP, V is orthogonal to and 
n^^. Hence, we can take V as the cross product 

V = n^"" X 

and W as a particular solution (e.g. the solution of mini- 
mum norm) of the overdetermined system of equations 

W"^n^ = 0, 
wT„SF^YTnSF. 

Step 2 : since HL belongs to FFP we have 

(W + fV)^nL-/ = 0, VfeM, 

hence n*^ verifies the overdetermined system of equa- 
tions 

W^n^=f, (3) 
\~^n^ = 0. (4) 

with the constraint Hn'^lP = 1. The computation of n'^ 
can be done by the following two steps: 

1 . compute W = V X W and V the minimum norm so- 
lution of system ([S]!-©, which gives a parametriza- 
tion 

= V + fW, 

of all its solutions, where f is an arbitrary real. 

2. determination of t such that ||n'"|p = 1: this is done 
by taking the solution t of the second degree equa- 
tion 

Wf 2 + 2W'^Vf + V'^V -1=0, 

such that n3 > 0. The tilt and swing angles are then 
obtained as 

0£ = — arcsin Mo 1 0l = — arcsin — - — . 

cos 01 

Step 3 : since SL is the intersection of LP and SFP, the 
coordinate vector U of a particular point U on SL is ob- 
tained as the minimum norm solution of the system 

V^n^ = 0, 

where we have used the fact that W G SFP. 

Step 4 : the translation of S can be computed such that 

U belongs to SP, i.e. 

(U-S)^n^=0. 

If we only act on the third coordinate of S and leave the 
two others unchanged, then can be computed as 

is = ^ . 

«3 




Fig. 7 Position of planes SFPi and SFP2 delimiting the depth of field 
region. 

Remark 1 When we consider a true camera lens, the nodal 
points H,H' and the front standard center L do not coincide. 
Hence, the tilt and swing rotations of the front standard mod- 
ify the actual position of the PSLP plane. In this case, we use 
the following iterative fixed point scheme: 

1 . The angles and Qi are initialized with starting values 

and el 

2. At iteration k, 

(a) the position of PSLP is computed considering 0^' and 
el 

(b) the resulting Hinge Line is computed, then the posi- 
tion of FFP and the new values 0^^' and 0^+^ are 
computed. 

Point 2 is repeated until convergence of 0^ and 9^ up to a 
given tolerance. Generally 3 iterations are sufficient to reach 
the machine precision. 

3.2 Characterization of the depth of field region 

As in the previous section, we consider that L and the nodal 
points H and H' coincide. Moreover, L will be the origin of 
the global coordinates system. 

We consider the configuration depicted in Figure Qwhere 
the sharp focus plane SFP, the lens plane LP and the sensor 
plane SP are tied by the Scheimpflug and the Hinge rule. 
The depth of field can be defined as follows: 

Definition 1 Let X be a 3D point and A its image through 
the lens. Let be the disk in LP of center L and diameter 
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f /N, where is called the f-number. Let K be the cone of 
base ^ and vertex A. The point X is said to lie in the depth of 
field region if the diameter of the intersection of SP and K is 
lower that c, the diameter of the so-called circle of confusion. 

The common values of c, which depend on the magnification 
from the sensor image to the final image and on its viewing 
conditions, lie typically between 0.2 mm and 0.01 mm. In 
the following the value of c is not a degree of freedom but a 
given input. 

If the ellipticity of extended images is neglected, the 
depth of field region can be shown to be equal to the un- 
bounded wedge delimited by SFPi and SFP2 intersecting at 
HL, where the corresponding sensor planes SPi and SP2 are 
tied to SFPi and SFP2 by the Scheimpflug rule. By mentally 
rotating SFP around HL, it is easy to see that SP is translated 
through and spans the region between SPi and SP2. The 
position of SPi and SP2, the f-number and the diameter 
of the circle of confusion c are related by the formula 



Nc _ pi -p2 
f Pi+Pi' 



(5) 



where pi, respectively p2, are the distances between the op- 
tical center L and SPi, respectively SP2, both measured or- 
thogonally to the optical plane. The distance p between SP 
and L can be shown to be equal to 



2piP2 

Pi +P2 



(6) 



the harmonic mean of pi and p2- This approximate defini- 
tion of the depth of field region has been proposed by vari- 
ous authors (see [8 . 1]) but when the ellipticity of images is 
taken into account a complete study can be found in JJ] . For 
sake of completeness, we give the justification of formulas 
(|5]l and (|6]l in Appendix lA.il In most practical situations the 
necessary angle between SP and LP is small (less that 10 
degrees), so that this approximation is correct. 



Remark! The analysis in Appendix I A. II shows that the ra- 
tio in equation (|5]l does not depend on the direction used 
for measuring the distance between SP, SPi, SP2 and L. The 
only condition, in order to take into account the case where 
SP and LP are parallel, is that this direction is not orthog- 
onal to n^. Hence, by taking the direction given by n^, we 
can obtain an equivalent formula to (|5]l. To compute the dis- 
tances, we need the coordinate vector of two points Ui and 
U2 on SPi and SP2 respectively. To this purpose we consider 
Step 3 of Algorithm [T] in section ITTl if W is the coordinate 
vector of any point W of HL, each vector U' can be obtained 
as a particular solution of the system 

U'^n^ = 0, 
U'"^n' = W"^n'. 



Since SPi can be defined as 

SPi = |xe]R^ (X-U')"^n^ = o|, 

and ||n^|| = 1, the signed distance from L to SPi is equal to 

d(L,SPi) =Ui"^n^. 
So the equivalent formula giving the ratio ^ is given by 



T 



dmSPi)-d{L,SP2) 



c/(L, SPi)+£/(L,SP2) 

(Ui-U2)^IlS 



(Ui + U2)^nS 



(7) 



The above considerations show that for a given orienta- 
tion of the rear standard given by n^, if the depth of field 
wedge is given, then the needed f-number, the tilt and swing 
angles of the front standard and the translation of the sen- 
sor plane, can be determined. The related question that will 
be addressed in the following is the question: given a set 
of points ^ — {X', . . . how can we minimize the f- 
number such that all points of ^ lie in the depth of field 
region? 



3.3 Depth of field optimization with respect to tilt angle 

We first study the depth of field optimization in two dimen- 
sions, because in this particular case all computations can be 
carried explicitly and a closed form expression is obtained, 
giving the f-number as a function of front standard tilt angle 
and of the slope of limiting planes. First, notice that has a 
natural upper bound, since ([5} implies that 



N < 



f 



3.3.1 Computation off-number with respect to tilt angle 
and limiting planes 

Without loss of generality, we consider that the sensor plane 
has the normal — (0,0, 1)^. Let us denote by Q the tilt 
angle of the front standard and consider that the swing angle 
^ is zero. The lens plane is given by 



LP = |x e M^ x"^!!^ o| , 



where 

= (0, - sine, cos 0)^, 

and any collinear vector to n'^ x is a direction vector of 
the Hinge Line. Hence we can take, independently of 0, 

V=(1,0,0)T 
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SP SP 




Using equation d?) the corresponding f-number is equal to 



SFP 



N{e,a) = 



h~t2 



h+ti 



(8) 



where we have used the notation a = {01,02). Finally we 
have, for ^ 

,„ , (fli — fl2)sin0 / f\ 

A^(0,a)=sign0:^-A-L_l^___(:^). (9) 



2cos0 — (ai +a2)sin0 \c 



Remork 3 When 0=0, then SFP does not intersect PSLP 
and the depth of field region is included between two parallel 
planes SFPi given by 

SFPi = {XeM-\X3=Zi}, /=1,2, 

where zi and Z2 depend on the f-number and on the position 
of SFP. One can show by using the thin lens equation and 
equation (|5]l that the corresponding f-number is equal to 



Fig. 8 The depth of field region when only a tilt angle 9 is used. 

and a parametric equation of HL is thus given by 

HL= {XeR-\ 3t eR, X = W(0) + fV}, 

where W(0) is the coordinate vector of a particular point 
W{9) on HL, obtained as the minimum norm solution of 

W(0)"^n^ =0, 
W(0)TnL=/. 

Straightforward computations show that for 0^0 



W(0) 



Consider, as depicted in Figure |8] the two sharp focus 
planes SFPi and SFP2 passing through HL, with normals 
ni = (0,-1, fli)"^ andn^ = (0,-1, a2)'^, 

SFPi = |xeM^ (X-W(0))^ni =o|, 

SFP2 = |xe]R^ (X-W(0))^n^=o|. 

The two corresponding sensor planes SPi are given by fol- 
lowing Steps 3-4 in Algorithm 1 by 




SFP, = {XeM-\X3=f,}, /=1,2, 



where 



/ 



Noizuzi) 



Zl ^ Z2 f 



(10) 



Oi sin — cos 



/= 1,2. 



3.3.2 Theoreticol results for the optimizotion problem 

Without loss of generality, we will consider a set of only 
three non-aligned points ^ = which have to 

be within the depth of field region with minimal f-number 
and we denote by X' , X^ , X^ their respective coordinate vec- 
tors. 

The corresponding optimization problem can be stated 
as follows: find 

(0*,a*)=arg min N{B,a), (11) 
GM 

fle^(0) 

where for a given the set ^(0) is defined by the inequal- 
ities 



(X'-W(0))V >0, 1,2,3, (12) 
(X'-W(0))^n^ <0, /= 1,2,3, (13) 



meaning that X^,X'^,X^ are respectively under SFPi and 
above SFP2, and by the inequalities 

-(X'-W(0))V+/<O, /= 1,2,3, (14) 

'1 v2 v3 



meaning are in front of FFP. 
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Remark 4 Points behind the focal plane cannot be in focus, 
and the constraints (fl4l i are just expressing this practical im- 
possibility. However, we have to notice that when one of 
these constraints is active, we can show that ai ~ cot 9 or 
fl2 = cot0 so that A^(a, 0) reaches its upper bound ^. More- 
over, we have to eliminate the degenerate case where the 
points X' are such that there exists two active constraints in 
(fT4l i: in this case, there exists a unique admissible pair (a, 9) 
and the problem has no practical interest. To this purpose, 
we can suppose that > f for i— 1,2,3. 

For 9^0 the gradient of A^(a, 0) is equal to 

/ COt0-fl2 

„N 2 Sign y / 
V(a,0)=: k\ -cot0 + ai 

^ ' ' (2cot0-(fli+a2))2 1 a,-a2 

V sin^e 

and cannot vanish since ai = 02 is not possible because it 
would mean thatX\X^,X^ are aligned. This implies that a 
lies on the boundary of si/ (9) and we have the following 
intuitive result (the proof is given in Appendix lA.2t : 

Proposition 1 Suppose that > f, for i = 1,2,3. Then 
when N{a,9) reaches its minimum, there exists i\, i2 with 
h 7^ h such that 

(X'l -W(0))^ni =0, 

(X'2-W(0))^n^ = O. 

Remark 5 The above result shows that at least two points 
touch the depth of field limiting planes SFPi and SFP2 when 
the f-number is minimal. In the following, we will show that 
the three pointsX^X^ andX^ are necessarily in contact with 
one of the limiting planes (the proof is given in Appendix 
lA3t: 



Proposition 2 Suppose that the vertices {^'}j=i...3 verify 
the condition 

Then N{a, 9) reaches its minimum when all vertices are in 
contact with the limiting planes. 

Remark 6 If is small, then N{9, a) in ^ can be approxi- 
mated by 

7 



.;V(0,a) = sign0(ai — fl2)sin0 



2c 



(16) 



and the proof of Proposition|2]is considerably simplified: the 
same result holds with the weaker condition 



iX'xX^'ll >0. 



(17) 



In fact, an approximate way of specifying the depth of field 
region using the hyperfocal distance, proposed in ||4l, leads 
to the same approximation of N{9,a), under the a priori 
hypothesis of small 9 and distant objects, i.e. 3> /. This 
remark is clarified in Appendix lA.4l 



We will illustrate the theoretical result by considering sets 
of 3 points. For a set with more than 3 vertices (but being 
equal to the vertices of the convex hull of ^), the determi- 
nation of the optimal solution is purely combinatorial, since 
it is enough to enumerate all admissible situations where two 
points are in contact with one plane, and a third one with the 
other The value 0=0 also has to be considered because it 
can be a critical value if the object has a vertical edge. We 
will also give an Example which violates condition (fTSl) and 
where A^(a, 0) reaches its minimum when only two vertices 
are in contact with the limiting planes. 

3.3.3 Numerical results 

In this section, we will consider the following function, de- 
fined for 9^0 



n{B) 



- min N(a,9). 

ae.c/ie) 



Finding the minimum of this function allows one to solve the 
original constrained optimization problem, but considering 
the results of the previous section, n(0) is non-differentiable. 
In fact, the values of for which n(0) is not differentiable 
correspond to the situations where 3 points are in contact 
with the limiting planes. We extend «(0) by continuity for 
= by defining 



«(0) = 



J i_ 

~"2 Zl 



_ J. _ 1 

Zl ^ Z2 f 



where 



zi = max Xl, Z2= min Xl . 

(=1,2,3 (=1,2,3 





This formula can be directly obtained by using conjugation 
formulas or by taking = in equation (l24l ). 

Example 1 We consider a lens with focal length/ = 5.1 O^^m 
and a confusion circle c = 3.10^^m (commonly used value 
for 24x36 cameras). The vertices have coordinates 



X^ 



Figure |9] shows the desired depth of field region (dashed 
zone) and the three candidates hinge lines corresponding to 
contacts EijVk'. 

- £'i2V3, obtained when the edge [X' is in contact with 
SFPi and vertex X^ with SFP2. 

- £'i3V2, obtained when the edge [X' is in contact with 
SFP2 and vertex X^ with SFPi. 

- £23^1, obtained when the edge is in contact with 
SFPi and vertex with SFP2. 
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Fig. 9 Enumeration of the 3 candidates configurations for Example 1 . 



Contact 


e 


n(9) 


E12V3 


0.0166674 


9.52 


EuVi 


0.025002 


7.16 


E23V1 


0.0500209 


28.28 



Table 1 Value of 6 for each possible optimal contact and correspond- 
ing f-number n(0) for Example 1. 



The associated values of and n{6) are given in Table [T] 
which shows that contact seems to give the minimum 

f-number. Condition (fTSl l is verified, since 



1, 



X X^^ 



||Xi-X2|| ' ||Xi-X3| 
11X2 xX3 1 



1.34, 



11X2 -X3 I 



1.48. 



Hence, the derivative of n{9) cannot vanish and the mini- 
mum f-number is necessarily reached for 9* = 0.025002. 

We can confirm this by considering the graph of n{9) 
depicted on Figure [TO] In the zone of interest, the function 
is almost piecewise affine. For each point in the interior 
of curved segments of the graph, the value of 9 is such that 
the contact of ^ with the limiting planes is of type ViVj. 
Clearly, the derivative of «(0) does not vanish. The possible 
optimal values of 9, corresponding to contacts of type EjjVk, 
are the abscissa of angular points of the graph, marked with 




-0.2 0.0 0.2 

tilt angle B 

(a) 



0.4 0.6 



■i 



£23Vl , 



V1V2 



V1V2 




V,V2 




0- 

00 



0.02 0.03 0.04 0.05 0.06 

tilt angle 

(b) 



Fig. 10 Example 1 graph of n{e) for 6 e [-0.6,0.6] (a) and 6 6 
[0, 0.06] (b). Labels give the different types of contact. 



Contact 


e 


n{e) 


E12V3 


0.185269 


29.49 


E23V1 


0.100419 


83.67 


E13V2 


0.235825 


47.04 



Table 2 Value of 9 for each possible optimal contact and correspond- 
ing f-number n{0) for Example 2. 



red dots. The graph confirms that the minimal value of n(0) 
is reached for contact £12^3. 

The minimal f-number isequalton(0*)=7.16.Bycom- 
parison, the f-number without tilt optimization is «(0) = 
28.74. This example highlights the important gain in terms 
of f-number reduction with the optimized tilt angle. 

Example 2 We consider the same lens and confusion cir- 
cle as in Example 1 (/ = 5.10"V c = 3.10"^^) but the 
vertices have coordinates 




The object is almost ten times smaller than the object of the 
previous example (it has the size of a small pen), but it is 
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Fig. 11 Example 2 graph of n(9) for 9 e [-0.6, 0.6] (a) and 9 e [0, 0.3] 
(b). Labels give the different types of contact. 



also ten times closer: this is a typical close up configuration. 
The values of and «(0) associated to contacts of type EjjVk 
are given in Table |2] which shows that contact ZiiaVi seems 
to give the minimum f-number We have 



2| 



|Xi-X2 , 

11X2 X X3 



0.16, 



X^ X X^^ 
XI -X3 



: 0.1775527, 



0.16, 



||X2-X3|, 

showing that condition ( fTSl ) is still verified, even if the val- 
ues are smaller than the values of Example 1. Hence, the 
derivative of n(0) cannot vanish and the minimum f-number 
is reached for 9* = 0.185269. 

We can confirm this by considering the graph of n{9) 
depicted on Figure [TT] As in Example 1 , the derivative of 
n{9) does not vanish and the graph confirms that the mini- 
mal value of n{9) is reached for contact EnVj. 

The minimal f-number is equal to n{9*) = 29.49. By 
comparison, the f-number without tilt optimization is n (0) = 
193.79. Such a large value gives an aperture of diameter 
0.26mm, almost equivalent to a pin hole ! Since the max- 
imum f-number of view camera lenses is never larger than 
64, the object cannot be in focus without using tilt. 



Fig. 12 Example 3 graphs of n(e) for 9 e [-1, 1]. 

This example also shows that Proposition|2]is still valid 
even if the object is close to the lens and the obtained opti- 
mal tilt angle 0.185269 (10.61 degrees) is large. 

Example 3 In this example we consider an extremely flat 
triangle which is almost aligned with the optical center From 
the photographer's point of view, this is an unrealistic case. 
We consider the same optical parameters as before and con- 
sider the following vertices 









:-•) 







for h 



0.01 and we have 
llXi xX2ii 



l|Xi-X2| 



= 0.02, 



11X2 xX^ 



|X> xX^I 
|Xi-X3| 



= 0.07, 



= 0.01, 



||X2-X3|, 

hence condition (fTSl l is violated in configurations V1V2 and 
ViVs- Since this condition is sufficient, the minimum value 
of n{9) could still occur for a contact of type EijVk- How- 
ever, we can see in Figure [T2K that n{9) has a minimum at a 
differentiable point in configuration V\Vj,. 
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Remark 7 Condition ( fTSb is not necessary: in the proof of 
Proposition [2] (given in Appendix IA.3l i. it can be seen that 
( fTSl l is a condition ensuring that the polynomial p{9) de- 
fined by equation jZST l has no root in [— ;r/2, ;r/2]. However, 
the relevant interval is smaller because X'' and X'^ do not 
stay in contact with the limiting planes for all values of 9 in 
[—n/2,n/2], and because all vertices must be in front of the 
focal plane. Anyway, it is always possible to construct ab- 
solutely unrealistic configurations. For example, when we 
consider the above vertices with h = 0.005, then no edge is 
in contact with the limiting planes and vertices X ' and X^ 
stay in contact with the limiting planes for all admissible 
values of 9. The corresponding graph of n{9) is given in 
Figure fT2b. 



3.4 Depth of field optimization with respect to tilt and 
swing angles 

As in the previous section, we consider the case of a thin 
lens and where the optical and sensor centers coincide re- 
spectively with the front and rear standard rotation centers. 
Without loss of generality, we consider that the sensor plane 
has the normal = (0,0, 1)^. The lens plane is given by 



LP = |xeM^ x"^ 



where 



n*" = (— sin^ cos 0, — sin 0, cos <j) cos 9)^ . 

A parametric equation of HL is given by 

HL= {XGM^ 3teR, X = W(0,0)+fV(0,0)}, 

where the direction vector is given by 

V(0,0) =n^ xn^ = (-sin0,sin0cos0,O)^, 

and W(0,0) is the coordinate vector of a particular point 
W{9,(j)) on HL, obtained as the minimum norm solution of 

W(0)^n^ =0, 
W(0)TnL=/. 

Consider, as depicted in Figure |7] the two planes of sharp 
focus SFPi and SFP2 intersecting at HL, with normals 
and respectively. The point W{9,^) belongs to SFPi and 
SFP2 and any direction vector of SFPi n SFP2 is collinear to 
V(0,0). Hence, we have 

SFPi = |xe]R^ (x-w(0,0))V =o|, 

SFP2 = |XGM^ (X-W(0,0))^n^ =0}, 



and 




Fig. 13 Position of the hinge line when both tilt and swing are used 
and = (0,0, 1)^. The LP and FFP planes have not been represented 
for reasons of readability. 



Using equation (|7]) the f-number is equal to 



N{9,0yy) = ^ 



(ui-u2)V 



(Ui+U2)^nS 



(19) 



where each coordinate vector U' of point U', for / = 1,2, is 
obtained as a particular solution of the following system: 

U'^n^ = 0, 

Remarks When — (0,0,1)^ the intersections of HL 
with the {L,Xi) and the {L,X2) axes can be determined, as 
depicted in Figure[T3] In this case, it is easy to show that the 
coordinates of the minimum norm W(0, 0) are given by 



W(0,(?)) = - 



/ 



sin cos^ 9 + sin 9 



(sin0cos0,sin0,O)^. 



We note that up to a rotation of axis (0,0, 1)^ and angle a, 
we recover a configuration where only a tilt angle y/ is used, 
where these two angles are respectively defined by 



sin y/ = sign 9 \J sin^ cos^ 9 + sin^ 9 , 



sma = 



(n^ xn^) X V(0,0) =0. 



sin ^ cos 9 
sini/ 



(18) 
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3.4.1 Optimization problem 

Consider a set 4 non coplanar points 2^ = 4, which 

have to be within the depth of field region with minimal f- 
number and denote by {X'},= i . .4 their respective coordinate 
vectors. The optimization problem can be stated as follows: 
find 

(0*,0*,ni*,n^*) = arg min A?(0, ^,n\n^), (20) 

e,<|i,nl,n2 

where the minimum is taken for , , and such that 
equation ( fTSl l is verified and such that the points i 4 

lie between SFPi and SFP2. This last set of constraints can 
be expressed in a way similar to equations (fT2t-(fT3]l. 

3.4.2 Analysis of configurations 
Consider the function n{9,(p) defined by 

n(0,(/)) = minA?(0,(/),n\ii^), (21) 

where and are constrained as in the previous section. 
Consider the two limiting planes SFPi and SFP2 with nor- 
mals and satisfying the minimum in equation ( |2TI ). Us- 
ing the rotation argument of Remark |8] we can easily show 
that SFPi and SFP2 are necessary in contact with at least 
two vertices. However, for each value of the pair {9,(j)), we 
have three types of possible contact between the tetrahedron 
formed by points {X'},=i. .4 and the limiting planes: vertex- 
vertex, edge-vertex, edge-edge or face-vertex. These config- 
urations can be analyzed as follows: let us consider a pair 
(00, 0o) and the corresponding type of contact: 

- Vertex-vertex: each limiting plane is in contact with 
only one vertex, respectively V, and Vj. In this case, n{0,^) 
is differentiable at (0o,0o) and there exists a curve 7,,,, 
defined by 

7„,(f)- (0(0,0(0), 7vv(0) = (00,00), 

such that ^n(0(f), ^(f)) exists and does not vanish for 
f 0. This can be proved using again the rotation argu- 
ment of Remark|8] Hence, n(0o, 0o) cannot be minimal. 

- Edge-vertex: one of the two planes is in contact with 
edge Eij and the other one is in contact with vertex Vk- 
In this case there is still a degree of freedom since the 
plane in contact with Eij can rotate around this edge in 
either directions while keeping the other plane in contact 
with Vk only. If the plane in contact with Ejj is SFPi, its 
normal can be parameterized by using a single scalar 
parameter f and we obtain a family of planes defined by 

SFPi(0 = |xeM^ iii(f)^(x-x') = o|. 



Contact 
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EnEu 


0.021430 


-0.028582 


12.35 


E21EU 


0.150568 


-0.203694 


90.75 


EnE2A 


0.030005 


-0.020010 


8.64 







0.100167 


88.18 


^243^1 


0.075070 


-0.050162 


42.91 


^134^2 


0.033340 


-0.033358 


9.64 


■^124^3 


0.018751 


-0.012503 


10.76 



Table 3 Value of 6 and (p for each possible optimal contact and corre- 
sponding f-number n{9,<p). 

For each value of f , the intersection of SFPi (f ) with PSLP 
defines a Hinge Line and thus a pair (0(f), 0(f)) of tilt 
and swing angles. Hence, there exists a parametric curve 

7,„(f) = (0(f), 0(f)), 7.v(0) = (0o,0o), (22) 

along which «(0,(/)) is differentiable. As we will see in 
the numerical results, the curve Yev{t) is almost a straight 
line when and (j) are small, and ^n(0,0(f)) does not 
vanish for f = 0. 

- Edge-edge: the limiting planes are respectively in con- 
tact with edges Eij, E^i connecting, respectively, vertices 
Vi,Vj and vertices V^jVi. There is no degree of freedom 
left since these edges cannot be parallel (otherwise all 
points would be coplanar). Hence, n{9,(p) is not differ- 
entiable at (00, 0o)- 

- Face-vertex: the limiting planes are respectively in con- 
tact with vertex V) and with the face Fiji; connecting ver- 
tices Vi,Vj,Vk. As in the previous case, there is no de- 
gree of freedom left and n{9,^) is not differentiable at 
(00, 0o)- 

We can already speculate that the first two configura- 
tions are necessary suboptimal. Consequently we just have 
to compute the f-number associated with each one of the 7 
possible configurations of type edge-edge of face-vertex. 

3.4.3 Numerical results 

We have considered the flat object of Example[T] translated 
in plane Xi ~ —0.5, and a complimentary point in order to 
form a tetrahedron. The vertices have coordinates 

All configurations of type edge-edge and face-vertex have 
been considered and the corresponding values of 0,0 and 
«(0.0) are given in Table [3] The £'i3£'24 contact seems to 
give the minimum f-number. Figure [14] gives the graph of 
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Fig. 14 Graph of «(e,(|)) for (9,(1)) e [-0.175,0.225] x 
[-0.3225,0.2775]. 















1 1 1 , , , , 




/ 




, • 




IV'iV, 

-/ / KRy / ^ ' 


1,1 




/ ,/ / iBss^iiL- 

- / / ^-^ITjWl. . 







-0.3 -0,2 -0,1 0,0 0,1 0,2 



n{9,(p) and its level curves in the vicinity of the minimum 
are depicted in Figures[T5]and[T6] In the interior of each dif- 
ferent shaded region, the (6,0) pair is such that the contact 
of ^ with the limiting planes is of type ViVj. The possi- 
ble optimal (0,0) pairs, corresponding to contacts of type 
EjjEiii of Fiji^Vi, are marked with red dots. Notice that the 
graph of n{9,(p) is almost polyhedral, i.e. in the interior of 
regions of type ViVj, the gradient is almost constant and does 
not vanish, as seen on the level curves. If confirms that the 
minimum cannot occur in these regions, as announced in 
Section 3.4.2. 

The frontiers between regions of type ViVj are curves 
corresponding to contacts of type EijVk and defined by Equa- 
tion (I22I 1. The extremities of these curves are (0,0) pairs 
corresponding to contacts of type EijEiii or Fiji^Vi. For ex- 
ample, in Figure [16] the (0,0) pairs on the curve separating 
V2V3 and V3V4 regions correspond to the £24 V3 contact. The 
extremities of this curve are the two (0,0) pairs correspond- 
ing to contacts -F124V3 and E\t,E2a- Along this curve, n(0, 0) 
is strictly monotone as shown by its level curves. 

Finally, the convergence of its level curves in Figure [T6] 
confirms that the minimum of n(0,0) is reached for the 
£'i3£'24 contact. Hence, the optimal angles are (0*,0*) = 
(0.030005, -0.020010) and the minimal f-number is equal 
to «(0*,0*) = 8.64. By comparison, the f-number without 
tilt and swing optimization is n(0,0) = 28.74. This exam- 
ple highlights again the important gain in terms of f-number 
reduction with the optimized tilt and swing angles. In our ex- 
perience, the optimal configuration for general polyhedrons 
can be of type edge-edge or face-vertex. 



Fig. 15 Level curves of n{9,<j>) and types of contact for {6,(j)) 6 
[-0.175,0.225] X [-0.3225,0.2775]. 















._-'^J'l24V3 












[135 


[EV:&lI 





^ — \ — -==^ ' 1 ' -■ — — ^ 1 — ' —I 1 

-0.035 -0.030 -0.025 -0.020 -0.015 -0.010 

Fig. 16 Level curves and types of of n(6,(j)) contact for {0,(j)) £ 
[-0.015,0.035] X [-0.0375,-0.0075]. 



4 Trends and conclusion 

In this paper, we have given the optimal solution of the most 
challenging issue in view camera photography: bring an ob- 
ject of arbitrary shape into focus and at the same time mini- 
mize the f-number This problem takes the form of a contin- 
uous optimization problem where the objective function (the 
f-number) and the constraints are non-linear with respect to 
the design variables. When the object is a convex polyhe- 
dron, we have shown that this optimization problem does 
not need to be solved by classical methods. Under realis- 
tic hypotheses, the optimal solution always occurs when the 
maximum number of constraints are saturated. Such a situa- 
tion corresponds to a small number of configurations (seven 
when the object is a tetrahedron). Hence, the exact solution 
is found by comparing the value of the f-number for each 
configuration. 

The linear algebra framework allowed us to efficiently 
implement the algorithms in a numerical computer algebra 
software. The camera software is able to interact with a robo- 
tised view camera prototype, which is actually used by our 
partner photographer. With the robotised camera, the time 
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Fig. 17 Construction of four approximate intersections of cones with 
SP 

elapsed in the focusing process is often orders of magnitude 
smaller than the systematic trial and error technique. 

The client/server architecture of the software allows us 
to rapidly develop new problem solvers by validating them 
first on a virtual camera before implementing them on the 
prototype. We are currently working on the fine calibration 
of some extrinsic parameters of the camera, in order to im- 
prove the precision of the acquisition of 3D points of the 
object. 

Acknowledgements This work has been partly funded by the Innova- 
tion and Technological Transfer Center of Region He de France. 



A Appendix 

A. 1 Computation of the depth of field region 

In order to explain the kind of approximation used, we have repre- 
sented in Figure[T7]the geometric construction of the image space lim- 
its corresponding to the depth of field region. Let us consider the cones 
whose base is the pupil and having an intersection with SP of diameter 
c. The image space limits are the locus of the vertex of such cones. 
The key point, suggested in 1 1 1 and 1 2 1, is the way the diameter of the 
intersection is measured. 

For a given line Si passing through L and a circle V of center L in 
LP let us call S) the set of cones with directrix S and base ^. 

For a given directrix S let us call A its intersection with SP, as depicted 
in Figure [TSli. Instead of considering the intersection of cones of direc- 
trix Si with SP, we consider their intersections with the plane passing 
through A and parallel to LP. By construction, all intersections are cir- 




SL 



(b) 

Fig. 18 (a) close-up of a particular intersection exhibiting the vertices 
of cones for a given directrix S. (b) geometrical construction allowing 
to derive the depth of field formula. 
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cles, and there exists only two cones Ki and K2 in ^C^, &) such that 
this intersection has a diameter equal to c, with their respective vertices 
A\,A2 on each side of SP, respectively marked in Figure [Tsl i by a red 
and a green spot. Moreover, for all cones in J(f{'ff, S>) only those with 
a vertex lying on the segment [^1,^2] have an "approximate" intersec- 
tion of diameter less that c. 

The classical laws of homothety show that for any directrix 3l, the 
locus of the vertices of cones Ki and K2 will be on two parallel planes 
located in front of and behind SP, as illustrated by a red and a green 
frame in Figure [T7t. Hence, the depth of field region in the object space 
is the reciprocal image of the region between parallel planes SPi and 
SP2 as depicted in Figure [T] 

Formulas JS) and (|6) are obtained by considering the directrix that 
is orthogonal to LP, as depicted in Figure 118b . If we note p = AL, 
pi = A\L, p2 = A2L, by considering similar triangles, we have 



Pi -P ^ P-P2 ^ 
Pi P2 / ' 



(23) 



which gives immediately 



P= ■ 



2piP2 



Pi +P2 

and by substituting p in {23), we obtain 

Pi -P2 ^ 

Pi +P2 / ' 

which allows to obtain (5). 
A. 2 Proof of Proposition [T] 



and straightforward computations give 



A'(a(e),e) ■- 



Xn J \ X^ X-i 



2cose- ( \ + \ )sine- 

, X^ X-, 



4+4 

Xn x^ 



(24) 

In order to prove the result, we just have to check that the derivative 
of N(a(e),e) with respect to 9 cannot vanish for 6 e [0, f ]. The total 
derivative of N{a{0), 9)) with respect to 9 is given by 



^«(a(e),9)^ 



■/(^T^)cos9+(^^-^jsin9^^ 



2 cos 9 - 



4 + 4 Isin9- 



4+4 

X, X3 



and its numerator is proportional to the trigonometrical polynomial 

p(9) =/?() + ^'1 cos 9 + ^2 sin 9, (25) 

where bo = X^' X'f - X^-X!^' ,b, = -f (x^^ - X^' ) , i<2 = /(^f - ^3' ) • 
It can be easily shown by using the Schwartz inequality that p(9) does 
not vanish provided that 



(26) 



Since Xp=x;' = 0, whe have &g = HXj' xX{' \\^ and bl + bl= f\\x{' 
Xj' Ip. Hence j26t is equivalent to condition )15t . this ends the proof. 



Without loss of generality, we consider that the optimal 9 is positive. 
Suppose now that only one constraint is active in l ll2t . Then there exists 
i'l such that (X'l — W) = and the first order optimality condition 
is verified: if we define 



g(a,9) = -(X' -W(9))^ni =X^' + 



/ 
sin 9 ' 



there exists Ai > such that the Kuhn and Tucker condition 

VA'{a,9) + AiVg(a, 9) =0, 
is verified. Hence, we have 



/ cot 9 — 02 \ 




-cot9+ai +Ai 1 


f -f \ 


I "l-fl2 / ' 


I co,s e / 


V sin- e 1 


V sin- 9 / 



and necessarily, a\ = cot 9 so that N[a, 9) reaches its upper bound and 
thus is not minimal. We obtain the same contradiction when only a 
constraint in l ll3t is active, or only two constraints in M2\ . or only two 
constraints in 1 113) . □ 



A. 3 Proof of Proposition |2] 

Without loss of generality we suppose that 9 > 0. Suppose that the 
minimum of ^(a, 9) is reached with only vertices i\ and !2 respectively 
in contact with liiniting planes SFPi and SFP2. The values of a\ and 
02 can be detennined as the following functions of 9 



A.4 Depth of field region approximation used by A. 
Merklinger 

In his book ((4|, Chapter 7) A. Merklinger has proposed the following 
approximation based on the assumption of distant objects and small tilt 
angles: if h is the distance from SFPi to SFP2, measured in a direction 
parallel to the sensor plane at distance z from the lens plane, as depicted 
in Figure[T2] we have 

2z H sin 9 

where H is the hyperfocal distance (for a definition see (5J p. 221), 
related to the f-number N by the formula 



H- 



Nc' 



and c is the diameter of the circle of confusion. Using i27\ . the f- 
number N can be approximated by 



N = sin 9 



'2zc' 



Since the slopes of SFPi and SFP2 are respectively given by ai and 02, 
we have 

h 

— — a\ — (325 

z 



X- 



and we obtain immediately 

N = [a\ — 02) sin 9 
which is the same as )16t . 
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Fig. 19 Depth of field region approximation using the hyperfocal. 
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